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<JJ ' Abstract 

We determine the radiative decay amplitudes for decay into D* and Z>y, or D* 
and Z? s 7 of some of the charmonium like states classified as X,Y,Z resonances, plus 
' some other hidden charm states which are dynamically generated from the interaction 

■ of vector mesons with charm. The mass distributions as a function of the -D7 or D s "f 

(f-) , invariant mass show a peculiar behavior as a consequence of the D*D* nature of these 

states. The experimental search of these magnitudes can shed light on the nature of 
CN ' these states. 



1 Introduction 



The use of the chiral unitary approach, combining chiral dynamics and unitarity in cou- 
pled channels has allowed one to study the interaction of pseudoscalar mesons, and of 
pseudoscalar mesons with baryons, at higher energies than allowed by perturbation the- 
ory [1]. One of the peculiar findings of the approach is that many resonances appear as 
poles in the scattering matrix as a consequence of the interaction, which are called dynam- 
ically generated resonances, and account for many of the low lying scalar meson and axial 
vector states, as well as for the low lying baryonic resonances. The success of this theory, 
providing properties of these states and accurate cross sections in production reactions, 
has stimulated the extension to the interaction of vector mesons. 

A natural extension of the chiral Lagrangians to incorporate vector mesons and their 
interaction is provided by the hidden local gauge formalism for vector interactions with 
pseudoscalar mesons, vectors and photons [2-5] , which provides a consistent and successful 
scheme to address many issues of hadron physics. Yet, as was the case with the interaction 
of pseudoscalar mesons or pseudoscalar mesons with baryons, also here it is the combina- 
tion of the interaction provided by these Lagrangians with unitary techniques in coupled 
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channels that allows one to obtain a realistic approach to the vector-vector interaction. In 
this direction, the work of [6, 7] has allowed to study the vector- vector interaction at inter- 
mediate energies, up to about 2000 MeV, where the nonperturbative unitary techniques 
are essential since many resonances are generated as a consequence of the interaction. In 
practice one solves a set of coupled channels Bethe Salpeter equations using as kernel the 
interaction provided by the hidden gauge Lagrangians and regularizing loops with a natural 
scale [8]. The results of [6] show that the / (1370) and / 2 (1270) mesons are dynamically 
generated from the pp interaction. Actually, there are strong experimental arguments to 
suggest that the /o(1370) is a pp molecule [9,10]. 

The work of [6] has been extended to the interaction of all members of the vector nonet 
with the result that eleven resonances are dynamically generated, most of which can be 
associated to known resonances, while other ones remain as predictions [7]. 

Some predictions of this approach for physical processes involving these states have 
readily followed to further support their nature as dynamically generated. In this sense the 
radiative decay of the / (1370) and / 2 (1270) mesons into 77 [11], were found in agreement 
with the experimental data. Similarly the J ftp decay into <p(uj) and one of the / 2 (1270), 
/ 2 (1525), /o(1710) resonances, and into K* and the 1^(1430) [12], was also found consistent 
with experiment. In the same line, the J/ip radiative decay into 7 and one of these 
nonstrange resonances was also found in agreement with experimental data [13]. More 
recently, the work of [11] has been extended in [14] to study the 77 and 7- vector meson 
decays of the eleven dynamically generated resonances of [7], with also good agreement 
with experiment in the cases that there are data. 

In [15] the extension to include charm mesons has been done studying the interaction 
of the p meson with D* mesons, where three states are obtained, one which can be easily 
associated with the tensor state D^AGO), another one which is very likely to be the 
L>*(2640) in view of its mass and the natural explanation for the small width compared 
with that from Z?2(2460), and a third one which corresponds to a scalar meson, for which 
no counterpart is yet reported in the PDG [16]. 

More recently the work has been extended to the interaction of D*D* in [17], where 
five resonances are dynamically generated, three of which could be tentatively associated 
to some X, Y, Z resonances reported recently, concretely the Y(3940), Z(3930), X(4160). 

Independently, an alternative approach to the hidden gauge formalism, based on chiral 
symmetry and heavy quark symmetry, but only in one channel, has been used in [18], 
where also bound states of the D*D* systems are found in some cases. 

Following the idea of [18, 19] that a Y(3930) and the Y(4140) in [17] could be actually 
D*D* and D*D* molecules, respectively, an idea was given in [20] that the shape of the 
spectrum in the radiative decay of these resonances into D*D^^ or _D*/) S 7, respectively, 
can further test the molecular assignment of these two resonances. 

We follow the idea of [20] with a different technical approach, and for the five states 
dynamically generated in [17]. The work of [17] provides scattering amplitudes for D*D* 
and its coupled channels. From there, evaluating the residues at the poles, one determines 
the coupling of the resonances to the different channels and this is all that one needs to 
evaluate the radiative decay in the Feynman diagrammatic approach that we follow. This 
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allows us to determine not only shapes of the spectra but also absolute numbers for the 
radiative decay in terms of the D*D^y coupling that can be taken from the experiment in 
the case of D* + — > D + ~f (and D*~ — > D~"y), and ratios of the R —>■ D*D^f decay width to 
the radiative width of the D*(D*) states in general. In [20] a different method based on 
wave functions of the states is reported and no absolute values are provided. In addition 
we give arguments on why the X — > D*D^ distribution with respect to the -D7 is the 
observable which connects easier with the dynamically generated nature of the resonances 
(molecular nature in the wave function picture). 

2 Formalism 

In [17] a coupled channel formalism was considered in which one had essentially the hidden 
charm D*D*, D*D* pairs plus all the charmless vector-vector pairs like pp, pu or 00, which 
have the same quantum numbers of the states that are investigated and provide the decay 
width of the XYZ states obtained. Five heavy states were generated, additional to the light 
ones obtained from the light vector pairs in [7], three of which were identified with states 
observed at Belle and Babar, the Y(3940), Z(3930) and X(4160), and two other states, so 
far not observed, which were called 3^,(3945) and 1^,(3912). The quantum numbers of the 
states and their assumed experimental counterparts are summarized in Table [U 



I G [J PC ] 


Theory 


Experiment 


0+[0++] 

0- [i+-] 

0+[2++] 
0+[2++] 

1- [2++] 


Mass [MeV] 
3943 

3945 
3922 
4157 
3912 


Width [MeV] 
17 



55 

102 

120 


Name 
Y(3940) 

'%(3945)" 
Z(3930) 
X(4160) 

'%(3912)" 


Mass [MeV] 
3943 ± 17 
39143i|j 

3929 ± 5 
4156 ±29 


Width [MeV] 
87 ±34 
33lf 

29 ± 10 

1391^ 


J PC 
J p+ 

2++ 



Table 1: Comparison of the mass, width and quantum numbers with the experiment. 

In [17] the states were identified by observing poles in the vector- vector scattering 
matrix with certain quantum numbers. The real part of the pole position provides the 
mass of the resonance and the imaginary part half its width. In addition the residues at 
the poles provide the product of the couplings of the resonance to the initial and final 
channels, from where, by looking at the scattering amplitudes in different channels, we can 
obtain the coupling of the resonance to all channels up to an irrelevant global sign, which 
is assigned to one particular coupling. In Tables [2] and [31 the couplings to the channels are 
also shown. As one can see from these tables, the states obtained correspond to basically 
bound D*D* or D*D* states, hence the decay into these pairs is forbidden, whereas the 
light vector-light vector channels provide the width of the states. However, if one looks at 



3 





Y(3940) 


y p (3945) 


Z(3930) 


X(4160) 


channel 




\m\ 


(MeV) 




D*D* 


18822 


18489 


21177 


1319 


D* S D* 


8645 


8763 


6990 


19717 


K*K* 


15 


40 


44 


87 


PP 


52 





84 


73 


ujijj 


1368 





2397 


2441 


# 


1011 





1999 


3130 




422 





1794 


2841 




1445 





3433 


2885 




910 





3062 


5778 


UJ<f) 


240 





789 


1828 



Table 2: Modules of the coupling in the I = 0, J = 0, 1, 2 sectors. 



the decay channel of the D* into D7, then the process X — > D*D^i is allowed, since the 
mass of the resonance X, for all the cases listed in Table dj exceeds the sum of masses of the 
final state. In Fig. [1] the corresponding Feynman diagram to the X — > D* + D~'~f process 
is shown. The D*~ propagates virtually between the production point of X — ► D* + D*~ 
and the decay point of D*~ — > D^'-f. This propagator is the relevant characteristic of 
the X — > D* + D~^ decay. Thus, this diagram is peculiar to the assumed nature of the 
resonance X as a molecule of D*D* and should be largely dominant over other possible 
processes [20]. The evaluation of this Feynman diagram is easy. All one needs is the 
coupling of the resonance to D* + D*~, together with the corresponding spin projection 
operator, and the vertex accounting for the decay of D*~ into D~^y. 

The spin projection operators on J = 0, 1, 2, evaluated assuming the three momenta of 
the D* and D* to be small with respect to the mass of the charmed vector mesons, which 
is indeed the case here, are given in terms of the polarization vectors by 

-p(O) = 1(1) (2) (3) (4) 

r 3 J * J ' 

V (l) = I(V(l) e (2) _ e W e (2hIf e (3) e (4) _ MM) 

x {i^+W-^S?^}- (1) 
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F p (3912) 


channel 


\9i\ (MeV) 


D*D* 


20869 


K*K* 


152 


PP 





puo 


3656 


pJ/iP 


6338 


p4> 


2731 



Table 3: Modules of the coupling in the I = 1, J = 2 sector. 



X{P) 




Figure 1: Decay of the X resonance to D* + D j. 

The amplitude obtained after summing all the diagrams included implicitly in the Bethe 
Salpeter equation, T = [1 — VG]~ 1 V, goes close to a pole, as g^g^is — s p ), where 
<?j(j) is the coupling of the resonance to the channel and are the spin projectors 
over spin k = 0, 1,2 of Eq. ([1]), see [6]. This final amplitude is depicted in the diagram 
of Fig. [2j In this way, if we take the case of the Y(3940), with J = 0, the first vertex in 
the diagram of Fig. [T]is -^e-^e- 2 ' 1 g D *£,* Fi, where Fi is the isospin factor needed to change 
from the isospin basis, where the couplings are evaluated in [17], to the charge basis. In 
the case of D* + D*~ , we have Fj = 4g. In what follows, we will call g the coupling of the 
resonance to the VV state in isospin basis. 

On the other hand the anomalous vertex for the D* decay into D7 is given by 

- itn*^ = -igpv^^p^e^D*) fc a e^ 7 ), (2) 

where p, k are the momenta of the D*~ and 7 respectively. This amplitude gives rise to a 
width 

r B ^5 7 = ^zalv^iMl - ml) 2 . (3) 
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Figure 2: Representation of the T matrix obtained from the Bethe Salpeter Equation 
in [15]. 



Unfortunately, only the value for the radiative decay of the D*~ — > D~"f and of its 
positive state partner are known. In this case we will be able to provide an absolute value 
for the radiative decay width of the XYZ resonances. In the other cases we will give the 
ratio of the radiative decay of the resonance to that of the D*. The value of gpy-y for the 
D*~ — > decay is given by 

gpVi = 1.53 x 10' 4 MeV-\ (4) 

which can be easily deduced using Eq. (J3j) from the experimental value of the width 
r = 1.54 KeV. 

Let us begin with the decay of the Y(3940). This state has isospin zero and spin zero. 
According to [17] it couples mostly to D*D*, has a smaller coupling to D*D* and very 
small coupling to pairs of charmless vectors, see Table [2J The couplings in [17] are given 
in isospin basis. However, we need them now in charge basis, which are readily obtained 
for the isospin combinations 

\D*D*,I = 0,I 3 = 0) = -^=(D* + D*- + D*°D*°), 

\D*D*, I = 1, J 3 = 0) = -1=(D* + D*- - D*°D*°), (5) 
\D s D s ,I = 0,I 3 = 0) = D*+D*~. 

Thus, the couplings of [17] for D*D* must be multiplied by l/y/2 to get the appropriate 
coupling for the charged or neutral states (a sign is irrelevant for the width), and do not 
require an extra factor for the case of D*D*. 

With the previous information we can already write the amplitude for the decay of the 
Y(3940) into D* + D~^, which is given by 

., (l) (2) i 

— it = —i—=g—=e\e 



V2 VS 1 1 P 2 -M 2 Dt +tM D *T D * 
XH)SPVW-/3^ 2 )A^( 7 ), (6) 

where the indices (1), (2) indicate the D* + and the D*~ respectively. The sum over the 
intermediate D*~ polarizations can be readily done as 

£4V< 2 > = -tf = -^, (7) 
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where we have neglected the three momenta of the intermediate D*~ which is in average 
very small compared with the D*~ mass, particularly at large invariant masses of the D~^( 
system which concentrates most of the strength, as we shall see. The sum of \t\ 2 over the 
final polarizations of the vector and the photon is readily done and, neglecting again terms 
of order p 2 /M|,*, we get the result 



£1*1 



1 1 

32 

1 1 

62 



~2 2 

9 Qpv-y 



~2 2 

9 9p Vl 



1 



P 



M 2 D , + iM D «T D * 



2 2 

p — m D 



p 2 — Mp* +iM D »T D * 



(8) 



The differential mass distribution with respect to the invariant mass of the D 7 system 

ly gi 

dT R 



M inv , with M 2 



p 2 , is finally given by 



1 



1 



dM inv AM 2 R (2tt) : 



■p*pu£I*I 



(9) 



where p* is the momentum of the D* + in the rest frame of the resonance X and po is the 
momentum of the D~ in the rest frame of the final D~7 system given by 



P 

Pd 



\'I\M 2 W M 2 D ^M\ 



2 ) 

inv 1 



2M A 
ML - m 2 D 



R 



(10) 



In the case of the tensor and spin one states we must do extra work since the projector 
operators are different. In this case we must keep the indices i,j in t and multiply with 
t* with the same indices i, j and then perform the sum over the indices This sums 
over all possible final polarizations but also the initial X polarizations, so in order to take 
the sum and average over final and initial polarizations, respectively, one must divide the 
results of the J2i,j tt* by (2 J + 1), where J is the spin of the resonance X. The explicit 
evaluation for the case of the tensor states, J = 2, of D*D* proceeds as follows: The t 



matrix is now written as 



-^9 9pv~i 1 2 \ e i e j +e i e i 



X 



p 2 -m 2 d , +iM D *r D , 



WP^ 2 W( 7 ). 



11^ 



As mentioned above, we must multiply tjj by t*j, recalling that the indices i,j are 
spatial indices and divide by (2 J+ 1) (5 in this case) in order to obtain the modulus squared 
of the transition matrix, summed and averaged over the final and initial polarizations. 
Neglecting again terms that go like p 2 /m*^ we obtain the same expression as in Eq. (jSJ). 
It is also easy to see that this is again the case for the J = 1 states. The normalization of 
the spin projection operators in Eq. makes this magnitude to be the same in all cases. 
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3 Convolution of the dY /dMi nv due to the width of the 
XYZ states 

Some of the dynamically generated XYZ states have a non negiglible width and, as a 
consequence, a mass distribution. That means there is a probability of these states to have 
a mass over the nominal mass and if one consider this fact, the PV7 decay width should 
increase. In order to consider this, we convolute the dT/dMi nv function over the mass 
distribution of the resonance R. We take Y/2 to both sides of the peak of the resonance 
distribution which account for a large fraction of the strength and produces distinct shapes 
in the 7D mass distribution. We find: 

1 r( M R+ T / 2 ) 2 ~ 1 1 

dT conv{r/2) /dM inv = — / dM 2 ( )Im-* dT/dM inv (12) 

1 mv N J(M R -r/2) 2 K tt j M 2 - M 2 R + iTM R ' K J 

with 

N = dM 2 ( — Jm- , 13 

J{m r -y/2Y n M 2 -M 2 R + iTM R 

As we will see in the next section, the use of Eq. (1121) leads to an increase of T(R — > PV'y) 
with respect to the result with the nominal mass M R . 



4 Results 

We show here the results for different cases: 



4.1 The Y(3940): Decay mode L>* + L> 7 

The results are the same reversing the signs of the charges. 

In Fig. [3] we show the distribution of Eq. (Q, together with Eq. ([8]), between the 
limits of M inv : mp and M R — mo*- Also, in order to see the effects produced when one 
considers the width of the state, we plot in the same figure dT conv /dM inv , taken from Eq. 
( |T2l) . We can see a very distinct picture, with most of the strength accumulated at the 
maximum values of M inv . The propagator of the intermediate D*~ and the factor (p.k) 2 
are responsible for that shape. In fact we show superposed in the same figure the result 
obtained ( normalized to the same area) substituting the propagator by a constant and 
removing the factor (p.k) 2 (or equivalently the factor (p 2 — m 2 D ) 2 ). We call Q the resulting 
distribution (or Q conv when one convolute this function taking into account the width of the 
R state). We can see that the pictures of dY/dM inv and Q (or equivalently dT conv /dM inv 
and Q conv ) are radically different and the reason is mostly due to the presence of the D*~ 
propagator which carries the memory that the resonance Y(3940) is assumed to be a D*D* 
molecule. The effects of considering the convolution are also visible in this picture. Now, 
dT conv /dM inv spreads beyond M R — mo*, and there is some probability for the state to 



S 




1870 1880 1890 1900 1910 1920 1930 1940 

M inv [MeV] 



Figure 3: The F(3940) -> D* + D~ T . Comparison of dT/dM inv and dT conv(r ^ /dM inv , 
and Q conv ( r / 2 ) as a function of the D~^y invariant mass. 



9 



decay into PV'-f up to Mi nv = Mr + Y/2 — mo*, where T is the width of the state. Also 
in this case, the difference between dY conv /dMi nv and <5 conv is clearly visible. 

For the case of decay into D*°D 0/ y the matrix element is formally the same except that 
now we do not know the experimental radiative decay width of the D*°. In this case we 
divide the mass distribution of the D*°D 0/ y decay by the width of the D*° —>■ D 0/ y and plot 
the magnitude 

i dr R 

r_D*^D7 dM inv 



with 

_ MB, - ml 

2M D * 

In Fig. H]we show the results of the dT r/ dM inv Y D * Dj distribution and also we compare 
with dTfl™ / dM inv T £>* D-y We can see that the enlarged range of the mass distribution 
between the limits Mi nv = Mr — mo* and Mr + Y/2 — mjj. is responsible for an increase 
in r(F(3940) -> D*°D ^). 

4.2 The y p (3945) 

This state has zero width, and here we show the difference between dY r/ dM inv and Q in 
the case of Yp(3945) — > D* + D~^ in Fig. [5] to see the effect of the inclusion of the D* 
propagator in Eq. (jSj). As one can see, the shapes can be clearly distinguished. Also, in 
Fig. E] we show the curves for dY R / dM inv Y D * Dl for the case of the neutral charm mesons 
in the final state. 

4.3 The Z(3930) 

This state has a larger width compared with the Y(3940) and 1^,(3945) states of 55 MeV, 
and for this reason the picture here is very different than in those cases when one takes 
into account this width. Thus, one can see a big difference between dY r/ dMi nv and 
dYft nv /dMi nv , Q and Q conv , as shown in Fig. [71 The relatively large width of the res- 
onance taken (55 MeV) is responsible for the different shapes compared to Fig. [7J a). 
Similar results are obtained for dY / dM inv Y D * Dl for decay into D*°D ^. 

4.4 The y p (3912) 

This case is very similar to that of the Z(3930). The shapes of dY r/ dM inv and Q are very 
different (also for dY c ^ nv /dM inv and Q conv ) as one can see in Fig. [SJ Now the width is 
considerably larger compared to that in the previous cases, since Y = 120 MeV. Similar 
results are obtained for the case of 1^,(3912) — > D*°D 0, -f. 



1 1 
26 



x 



~2 2 

9 9py, 



p- 



m 



D 



48vr Ml, 



p 2 -M 2 D *+iM D *Y D » 
1 1 _ 



k{M 2 D , 



m 2 D ) 2 AM 2 (2tt) 3 



P Pd, 



(14) 



10 



0.00012 



> 0.00010 
4^-00008 

ufo.00006 

> 

jf 0.00004 
*3 0.00002 



0.00000 




1880 



1900 1920 
M inv [MeV] 



1940 



Figure 4: The F(3940) -> D*°D° r . Comparison of dT / dM inv T D * Dl and 

d r c ° nv ( r / 2 )/dM im) r D * D7 . 



n 




12 




Figure 6: The F p (3945) -> D*°D r . dT / dM inv T D * Dl and Q as a function of Mj. 
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Figure 7: The Z(3930) — > D* + D 7: a) dT/dM inv and Q as a function of M inv . b) 
dT conv /dM mv and Q conv . 

4.5 The X(4160) 

In this case the isospin factor is Fj = 1 rather than 1 / y/2. The formula is the same as 
before removing a factor 1/2 in Eq. ([8]). Once again we do not have the experimental 
decay rate for the radiative decay of D*~ and we plot the results for Eq. A14H in Fig. [91 In 
this case the decay into D* + D~ / y is also possible. However, the coupling to D*D* of this 
resonance (also assumed to be a D* + D*~ molecule in [20] and [19]) is found small in [17], of 
the order of 17 times smaller, hence the rate for this channel should be drastically smaller. 
In order to test the D* + D*~ component of this molecule, the allowed strong decay into 
D*D* is preferable. This latter measurement is a more efficient tool to get the strength of 
this coupling and compare with the theoretical predictions. 

In Tables H] and [5] we show integrated values for T(R — > PV'y) and also rates of T(R — > 
PVj) with respect to T(DJs — > D^j). In the case of the decays of the resonance into 
D*°D°~f, D* + D~j, D*°D 7i° and D* + DJtt°, we compute g PVl in Eq. © taking T(D*°) < 
2.1 MeV and T(D* + ) < 1.9 MeV. We show in Table H the integrated values for T(R -> 
D* + D~ / ~f) which are very small, of the order of 1CT 1 — 1 eV if one does not consider the 
convolution of the dT/dM inv distribution. However, when one considers the width of the 
XYZ resonances given in Tabled], these values become bigger (about one order of magnitude 
in some cases). 

In the case of the X(4160) we can only put a boundary for the T(X — > Dl + Dj , y), which 
is 39.9 KeV, but we give rates of T(X -> D* s + D~^) respect to T(D* S ~ -> D77) in Table 
HI For this observable we get a value of 2.3 x 10~ 2 and 0.14 before and after convolution 
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conv (r/2) J* 
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M lnv [MeV] 



\ 

\ 

N 

1940 1960 



Figure 8: The F p (3912) -> D* + D~r- a) dT/dM inv and Q as a function of M inv . b) 
rfr conv /rfM int , and Q conv . 
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Figure 9: The X(4160) -> D* + D~r- a) dT R / dM inv T D » Dsl and Q as a function of M inv . b) 
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State Decay 


T [keV] 




r conv(r/2) [ ke y] 


pconv(I72) /p 


Y(3940) D*+D-j 


2.7 x 10- 3 


1.8 x IO- 3 


2.9 x 10- 3 


1.9 x 10- 3 


Y p (3945) D*+D~^ 


3.1 x 1(T 3 


2.0 x 10~ 3 






Z(3930) D* + D-j 


4.1 x 10^ 4 


2.6 x 10~ 4 


1.0 x 10~ 3 


6.7 x 10" 4 


Y p (3912) D*+D~-f 


1.0 x 1(T 4 


6.7 x 10~ 5 


2.7 x 10~ 3 


1.8 x 10~ 3 


X(4160) D*+D 7l 


< 39.9 


2.3 x 10- 2 


< 2.4 x 10 2 


0.14 


Table 4: Decay of the XYZ resonances 


into D* + D-j and D* + D;j. 


State Decay 


T [keV] 


r/r D * _> DOl r conv ( r / 2 ) [keV] r 


conv(r/2) /r 


Y(3940) D*°D°7 


< 2.6 


3.2 x 10~ 3 


< 2.7 


3.4 x 10~ 3 


Y p (3945) _D*°D°7 


< 2.9 


3.6 x 10~ 3 






Z(3930) D*°D°-f 


< 0.48 


6.0 x 10~ 4 


< 1.0 


1.3 x 10~ 3 


F p (3912) D*°D°7 


< 0.15 


1.9 x 10~ 4 


< 2.4 


3.0 x lO" 3 



Table 5: Decay of the XYZ resonances into D*°D 0, -f. 

respectively. When the final state contains neutral charm mesons, we give both amplitudes 
and rates which can be seen in Table [5j In Table El we see that T/T £>.o^^o 7 is of the order 
of 10 4 — 10~ 3 for all the states before the convoluting dY r/ dM inv and becomes larger when 
one convolutes this function. 



5 Summary 

We have presented results for decay of the heavy dynamically generated states from the 
vector-vector interaction, with hidden charm, into D* and -D7, or D* and l) s 7- We find a 
very distinctive shape in the D7 and D s ^ invariant mass distributions, which is peculiar 
to the molecular nature of these states as basically bound states of two charmed vector 
mesons. It was suggested in [17] that some of these states correspond to some of the 
X,Y,Z states found at the Belle and Babar facilities. We hope the findings of the present 
paper stimulate experimental work in this direction to further learn about the nature of 
the X,Y,Z resonances. 
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